Electronic and optical trends in carbon nanotubes under pure bending 
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The high aspect ratio of carbon nanotubes makes them prone to bending. To know how bend- 
ing affects the tubes is therefore crucial for tube identification and for electrical component design. 
Very few studies, however, have investigated tubes under small bending well below the buckling 
limit, because of technical problems due to broken translational symmetry. In this Letter a cost- 
effective and exact modeling of singe-walled nanotubes under such small bending is enabled by 
revised periodic boundary conditions, combined with density-functional tight-binding. The result- 
ing, bending-induced changes in electronic and optical properties fall in clear chirality-dependent 
trend families. While the correct trends require full structural relaxation, they can be understood 
by one general argument. To know these trends fills a fundamental gap in our understanding of the 
properties of carbon nanotubes. 
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Carbon nanotubes (CNTs), with an aspect ratio akin 
to tens of meters long human hair, bend frequently. [U |2] 
In most transmission electron micrographs CNTs appear 
as curly web of hair, unless placed on a support with 
special care.|3tt6^ Bending is hard to avoid in experiments 
completely. 

With theory the state of affairs is just the opposite. 
On one hand, it's easy to model infinitely long and per- 
fectly straight tubes, using Bloch's theorem and periodic 
boundary conditions. On the other hand, it's usually 
difficult to model large-scale distortions once the trans- 
lational symmetry gets broken. Bending has been there- 
fore commonly modeled by finite tubes and thousands of 
atoms, accessible only by classical potentials. [7, 8 

However, since CNTs serve mostly as electronic com- 
ponents, bending and other distortion studies ought 
to incorporate the electronic structure. [9 Sure enough, 
density- functional theory has been used to study CNT 
bending, but reliable electronic trends are hard to ex- 
tract, because tubes are short and bending violent — 
mostly near or above the buckling limit. [10, 11 Results 
also depend on how the tubes are bent, whether bending 
is achieved by forces from external bodies, or by geomet- 
rical constraints. 

In this Letter a cost-effective way to model tubes un- 
der pure bending is enabled by revised periodic bound- 
ary conditions (RPBC). In this approach the bending 
emerges naturally, being caused by boundary conditions 
alone, and is devoid of spurious effects. While similar 
approach for carbon nanotubes was used earlier by Du- 
mitrica et a/.,[T2j[T3] as well as Malola et a/., [14], they did 
not study electronic properties. The focus in this work 
is, therefore, to answer the question: "How do band gaps 
and optical properties change as single-walled CNTs get 
slightly bent?" The answers to this question, as it will 
turn out, have clear trends that can be understood with 
one fundamental argument. 

To modeling of bending is done by replacing the usual 
translation symmetry by a rotation symmetry around a 
given origin, as clarified in Fig. [l^; for details of this 
RPBC approach, see Ref. [HI Since the rotation is a 
symmetry operation, the modeled system as a whole is, 
in fact, a huge nanotorus. The amount of bending is 



measured by the parameter 



where D is diameter and R torus's radius. This param- 
eter is universal, and enables comparing tubes with dif- 
ferent diameter. (The sidewall strain on the compressive 
side is e ^ —& and on the tensile side e ^ independent 
of D) I concentrate here on bending around the experi- 
mentally relevant O~0.01orl% region, with 5 % at 
maximum; [16 1 % corresponds to bending a human hair 
into a loop with diameter of ~ 1 cm. 

The electronic structure itself is modeled by density- 
functional tight-binding (DFTB) method, [17 using the 
hotbit code with an RPBC implementation. \T8\ fT9] This 
method uses first principles theory first to parametrize 
matrix elements, then to fit pair potentials that give to- 
tal energies. [18 DFTB describes carbon materials rea- 
sonably well, regarding both mechanical and electronic 
properties. [2QH22] Considering these literature records, 
DFTB should describe the physics of bend CNTs well. 
Note that the approach to treat pure bending is exact, the 
sole approximation in this work being DFTB itself. [23 

The simulation in practice ran as follows. I selected 
CNTs with different chiralities (n, m) such that the num- 
ber of atoms in the unit cell was below 210, which made 
some 60 tubes in total. All the tubes were then bent from 

% to 5 % in steps of 0.25 % in the following way: First, 

1 fixed a according to the O of interest, using reasonable 
estimates for the initial geometry (including estimate for 
R). Second, I optimized the tubes to the maximum force 
of 10"^ eV/A.[l4l[25] Note that the only parameter fixed 
was a; the tube can freely move in radial direction, to 
find the optimum geometry and radius — the bending is 
pure. The number of equally spaced /^-points (/c-point 
equivalent for rotational symmetry) were 50 A/L for op- 
timizations and 500 A/L for electronic structure analysis, 
where L is the unit length of the straight tube. While 
the simulation cells contained ~ 100 atoms, the largest 
corresponding torii would have contained ~ 10^ atoms, 
giving an obvious motivation for this approach. 

Prior to presenting the results for bent CNTs' elec- 
tronic structure, it's illustrating first to review bent 
tubes' structural trends. [26 Namely, the inset in Figflb 
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FIG. 1: (color online) (a) (10, 5) CNT, with D = 10.3 A, bent 
to G = 1 % and G = 5 %. Red (light gray) atoms denotes the 
simulation cell, and the symmetry operation is a rotation of 
an angle a around given origin. G is a universal measure of 
bending: the impression for the amount of bending for given 
G is independent of D. (b) Elastic bending energy (per unit 
length and divided by tube diameter) as a universal function; 
the dashed line is the analytical expression, Eq.(|3]). Right 
inset: cross-section view of atoms' positions from all tubes, 
plotted with DG = 0.3 A. The dashed line is an ellipse with 
r = 0.949, as given by Eq.H with D(c> = 0.3 A. The flatten- 
ing is only somewhat visible when comparing to the perfect 
square. Left inset: ratio of the two major axes of the flattened 
cross-section. Dashed line is the analytical expression 



shows that bent tubes get slightly flat. Flattening is due 
to two competing effects: flatter tubes have on one hand 
less energy due to strain, but on the other hand more 
energy due to increased sidewall curvature. Using thin 
sheet elasticity theory to calculate the optimum between 
these competing effects, one gets 



r = 1 



4 [(A/6L))2 - 1] 



(2) 



for the ratio of the minor and major axes of the resulting 
elliptical cross-section, where A = ^Jl^t^/Y ^ 1.16 A, 
K = 1.6 eV is graphene's bending modulus and Y ^ 
340 PaA its Young's modulus in two dimensions. Equa- 
tion Q shows that thicker tubes flatten more easily, and 
that tubes with equal DQ are equally flattened; the inset 
of Fig.[l]3 shows all tubes' cross-section with DQ = 0.3 A, 
or r = 0.949. Since buckling means complete flattening 
of the tube, there ought to be some critical limit Tc below 
which tubes buckle. Earlier calculations have shown that 
Tc ~ 0.65, which means DGc 0.75 A; [8, 26^ the tubes 
here have Oc = 5 . . . 20 %, and hence always Q < Qc 
The analytical expression for the elastic energy of tubes 



flattened this optimum way is 

E = ZdYQ^ 



(3) 



per unit length. [27^ The quantity E/D^ which hence is 
universal for all tubes, is plotted in Fig. [TJd. Small- 
diameter tubes deviate most from this behavior because 
thin-shell theory is less valid. This structural behavior of 
CNTs under bending agrees with previous classical simu- 
lations (using thousands of atoms) . [7^. [HI [TTl [26 and thus, 
for its part, validates the modeling approach. 

Now, let's finally turn the attention to the results in 
electronic structure. Figure |2^ shows the density of states 
and band structures for straight and bent (10, 0) tubes. 
The first observation is that the bending-induced changes 
in the electronic structure are gradual. Band extrema — 
the van Hove singularities — move gradually up and down, 
while some band anticrossings become more visible. The 
band gap increases little, yet distinctly. Figure [2b shows 
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FIG. 2: (color online) (a) and (b) Density-of states (left) and 
band-structure (right) for straight and bent (10, 0) and (11,0) 
tubes, (c) The behavior of band gap A(G) as a function of 
bending depends on tube's ^-family [defined in Eq.Q]. Gaps 
are renormalized by the gaps of the straight tubes, A|G = 0). 
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the same results for (11,0) tube. The changes due to 
bending appear similar, but note that the trends are just 
the opposite: the band extrema move in the opposite 
direction, and band gap decreases. 

The trends in the band gaps A(0) upon bending, sum- 
marized in Fig. [2]3, is the first main result of this work. 
The trends fall in families depending on CNT chirality 
(n,m) as 



a) 



q = {n — m) mod 3; 



(4) 



hence gap decreases for g = and q = 2^ whereas it 
increases for g = 1. For q ^ semiconducting tubes 
having sizeable gaps the relative change in A is smaller 
than for q = -tubes with a minigap {q = tubes are 
often referred to as being "metallic"). The q = tubes 
are therefore affected the most by bending — the bending 
can almost make those minigaps vanish. 

Trends above have two exceptions, armchair and small 
diameter tubes, that are excluded from Fig. First, 
the electronic structure of armchair tubes (n = m), it 
turned out, are very robust against bending. They re- 
mained metallic, with band structure nearly intact. This 
robustness, explained by symmetry arguments, has been 
reported earlier. [28, 29 Second, the electronic structure 
of tubes with small diameter (< 5 A) can change drasti- 
cally. Those changes are interesting, but they are unsys- 
tematic and impossible to understand jointly, and there- 
fore outside the scope of this paper. 

The next trends I shall discuss are optical trends. They 
are investigated via the imaginary part part of the di- 
electric function, £2(00)^ which is directly proportional to 
the optical absorption. The function 62(00) is calculated 
within the random phase approximation. [30^ This calcu- 
lation includes, however, one issue that I wish to discuss 
first. Assuming both incident and emitted light polar- 
ization in the Cartesian direction /3, the calculation will 
involve matrix elements {ipi^^a \p^\i^K,b)^ where p is the mo- 
mentum operator and where the states ip^^i are delocal- 
ized over the whole torus. [30 While the optical response 
could indeed be calculated for the nanotorus, the em- 
phasis here is on the optical properties arising from bent 
fragments of CNTs, not on nanotori as whole, complete 
objects. I therefore evaluate the matrix element only 
within the chemical interaction range from the unit cell; 
since the tube axis reorients very little between neighbor- 
ing unit cells, fixing the orientation of coordinates makes 
sense. The resulting 82(00) is hence a local property: the 
situation resembles probing and measuring the tube us- 
ing a laser with a spot size much smaller than the radius 
of curvature, as sketched in Fig. [3^. 

The function 82(00) calculated this way for a straight 
(10, 0) tube is then shown in Fig. [3|i, on the left. On the 
right, figure shows further how tne first two transition 
peaks Ell (between the first pair of occupied and un- 
occupied van Hove singularities) and E22 (between the 
second pair of occupied and unoccupied van Hove sin- 
gularities) change as a function of 9: £^11 increases and 
£^22 decreases in energy, as can be deduced also from the 
densities of states in Fig. [2^. Figure ^jp shows the same 
for (11,0) tube, and the behavior is, again, the opposite: 
£^11 decreases and E22 increases in energy. 
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FIG. 3: (color online) (a) Left: imaginary part of the dielec- 
tric function for a straight (10, 0) tube; incident and reflected 
light polarizations are either parallel (zz) or perpendicular 
(xx) to tube axis. Right: the effect of bending on the optical 
£11 and £22 peaks, (b) The same for (11,0) tube, (c) The 
behavior of £11 (left) and £22 (right) transitions depends on 
the tube ^-family, Eq.Q. Transition energies are renormal- 
ized by straight tube energies. The irregularities with q — 2 
are not real, but related to numerical problems in structure 
optimization; they indicate, however, the sensitivity to prop- 
erly optimized geometries. 



The trends in transitions £11 and £22 upon bending, 
summarized in Fig. (Sj), is the second main result of this 
Letter. Since A is essentially the same as £11, the left 
of Fig. [s]^ shows the same trends as Fig. (here £11 for 
(7 = is across the first transition in the optical range, 
not across the miniband). For £22, however, the trend is 
reversed: £22 increases ioi q = 2 and decreases ioi q = ^ 
and g = 1. For the trends in the absolute transition en- 
ergies, with different levels of tight-binding, see Refs. 1211 
and 1221 

The analysis of 82(00) excluded again armchair and 
small diameter tubes, armchair tubes for being insen- 
sitive for bending, and small diameter tubes for show- 
ing too unsystematic behavior. I further mention that 
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properties with transverse-polarized light (dotted lines 
in FigsjS^i and ^) are insensitive to bending; the tube 
cross-section remained, after all, rather spherical. 

At this point it's pertinent to mention that these trends 
appear to contradict some earlier works. Reference ST! re- 
ported bend-induced gap widening for g = and Refs. 29 
and I28l concluded that bending has no effect because — 
so the argument ran — opposite strains on opposite sides 
of the tube circumference cancel out. The very title of 
Ref. 32 suggested bend-induced insulating gap for zigzag 
and chiral CNTs. While my results do agree on the insen- 
sitivity of armchair tubes towards bending, [28, 29 why 
are there contradictions in those other trends? 

This question brings me to the third main result of 
this work: the correct trends are obtained only after full 
structural relaxation. The absence of relaxation in previ- 
ous works, therefore, explains the above contradictions. 
Indeed, if the calculations were repeated with a mere 
skilled estimate for the bent geometry, the trends re- 
garding (^-families turned out wrong. Even single-orbital 
tight-binding explains the trends correctly — if it's only 
given the properly optimized geometries. 

What makes the relaxation so important, is that, 
due to anharmonic effects, the carbon- carbon bonds will, 
on average, stretch upon bending. Pure stretching of 
CNTs, in turn, has been reported to yield the very same 
trends. [28^ In the limit 0^0 bonds do behave harmon- 
ically and opposite strains on opposite sides of the tube 
do have a cancelling effect; note that all curves in Figs|2]3 
and [Sj) start with a zero slope, that is, anharmonicity 
emerges only in the second order in 0. As bending in- 
creases, the bonds in the outer part stretch progressively 
more than bonds in the inner part shrink. The simulated 
circumference-averaged strains increase as Savg ^ 0.4-0^, 
yielding some 10~^ strains with = 5% bending. 

To conclude, the trends above ought to be correct 
and independent of any approximations of tight-binding, 
since the origin of the trends, the anharmonicity of 
carbon-carbon bonds, is so plausible and fundamen- 
tal. Measuring these trends from individual tubes is 
feasible, [6 although care should be taken to avoid other 
distortions and to make the bending pure. More rele- 
vant, however, are still the general trends, not properties 
of individual tubes — because we still don't have full ex- 
perimental control over nanotube chiralities. 
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